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A generalized non-local nonlinear Schrödinger equation describing the phenom- 
enon of nonlinear Landau damping in a relativistic two-component plasma is 
deduced using the kinetic-theory approach of Vlasov. Parameters appearing in 
the equation are evaluated explicitly for the case of a Maxwellian distribution 
function. 


1. Introduction 


The phenomenon of Landau damping occupies a central position in the study 
of many important properties of plasmas (Davidson 1992). In general, there are 
two kinds of Landau damping: linear and nonlinear (Weiland, Ichikawa & 
Wilhelmsson 1977; Weiland 1981; Mjolhus & Wyller 1986). Perhaps the first 
attempt to incorporate the effects of Landau damping in the study of solitons 
was that of Ott & Sudan (1969, 1970), who derived a nonlinear dispersion 
relation containing a term representing Landau damping. Several attempts 
have since been made to study the influence of this important phenomenon in 
other situations, such as in relativistic collisionless plasmas (Roy Chowdhury, 
Pakira & Paul 1988), magnetized plasmas (Som et al. 1992) and two-temperature 
systems (Chen (1984). Landau damping has been discussed in detail by Chen 
(1984). 

Ichikawa, Imamura & Taniuti (1972) studied the effect of resonant particles 
having group velocity v, = 0w/0K on finite-amplitude plasma waves. They 
showed that such waves can be described by a nonlinear Schródinger equation 
with an extra nonlinear non-local term that actually represents nonlinear 
Landau damping. A Lagrangian formalism for nonlinear wave packets was 
adopted by Dewar (1972), who determined a non-local nonlinear frequency shift 
and obtained a similar equation. All of these considerations were non- 
relativistic. Therefore here we consider a collisionless isotropic plasma consisting 
of electrons and ions, with the electrons being assumed to be relativistic. The 
whole system is described in the framework of Vlasov kinetie theory 
(Kadomtsev 1965; Baryakhtar et al. 1992). We have then followed the approach 
of Fried & Ichikawa (1973) to investigate the nonlinear dispersion relation, and 
hence have obtained a nonlinear Schróinger equation with an extra nonlinear 
non-local term. As mentioned above, this term represents nonlinear Landau 
damping in the relativistic plasma. 
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2. Formulation 


Let f, and f, denote the velocity distribution functions for the electrons and 
ions in the plasma, with the electrons assumed to be relativistic. With ¢ 
denoting the electrostatic potential, the Vlasov kinetic equations describing a 
collisionless two-component magnetized plasma are 


of, |, Se 9b fe _ 
etas дй ^ 


ôx т, 0x dv, = 


ZE = ane [(,—/)4 


We denote the mass of the electron by m; = m, (1—v2/c?)7}, m, being the rest 
mass. For the ions m; = m, since they are taken to be non-relativistic. 
Converting all the momentum variables to the corresponding velocities, we can 
rewrite the above set of equations as 


of... Ofe 2 d m 
t +v * Ox AT x дэл = 0, 
"RP 
at Se т Ox ди, | ý A 
ad _ 
aoc 476" f; dv, + 4nem, |f, nwe) dv,, 
2 
where qv 0) = 1s. 
_ 1—v2/2c? 
Хх) = FETTE 


To derive the dispersion relation, we set 
fe = Ew) + fei 
Л = Fo) + fa; (3) 
ф = hy, 
where F, and F, are the equilibrium values of the corresponding distributions, 


and fas / and ф, are perturbations, each of which is proportional to 
exp («KX — 101). Substitution of (3) into (2) leads to 


_ _ tato) Е) 

Ла mE mv, —w/K) ) d, 
_ _ Fw) 

fa = m,(v;— w/K) OA 
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v) dv. (Ve) Xv.) р, 
d 25 і 1 2 e e 
an 1 = 4пе? | em 6 е Kv, оу) 9) dv, (5) 
which will be written as 
e(K,w,0) = 


and is the linear dispersion relation not involving the amplitude. A method of 
computing corrections to (5) was given by Fried & Ichikawa (1973), who 
showed how a nonlinear Schrödinger equation can be deduced from it. In the 
following section we shall briefly discuss this method. 


3. Derivation of the NLS equation 


The basic idea behind the analysis of Fried & Ichikawa (1973) is to consider 
the form of the wave at time t = 0 to be 


$(Х,0) = [axe eX + С.С., (6) 


where фк is peaked around К = К, (С.С. stands for complex conjugate). In 
linear theory the long-time behaviour of the system is given by 


$(x,t) = fax фк ETD CC, (7) 


where w = w(K) is the least-damped solution of the dispersion relation. 

Here we shall consider the case where both the amplitude and the spread of 
K values are small, so that an expansion of the dispersion relation, first in A? 
and then in К — K,, is appropriate; thus 


(К, 9, A) = є(К, о, 0) 41 06+... =0. (8) | 
F ge 
Or о, we get w = O(K)-A EUM (9) 


where Q(K) is a solution of the linearized dispersion relation 
ЕК, Q(K), 0] = 0. (10) 
Again, Q(K) = QC K - Ky) = ov Ev, H+ 2H*+..., (11) 
with w, = Q(K,), v, = Q'(K,), v, = Q'(K,) and Ж = K —K,, where the prime 
denotes the derivative with respect to K. Thus we get 
w = Qo v | 
v= VA +w ЖА, | 


= So eas 
Ё 0433 24| ЕРУ 


(12) 


with 
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Substituting (6) into (7), we obtain 


$(z,1) = | IK фк, КО THOM О.С, 
= уба, t) екоо + C.C., (13) 
yz, t) = fax Parr i 2". (14) 

From (12) and (14), we can immediately infer that у satisfies 
(es E az — ру = о, (15) 
which is the usual nonlinear Schrödinger equation, when we have identified || 


as the amplitude A. However, the dependence of є on A? may be of a more 
general form 


[ew S(x — x') A(x’), 
where S(x) is a kernel, whence 


e(K,wA) = в(К, 0,0)+ [dx Sisi afe 


Now, дє/дА? = [ax вее), 
-1 
so that и = - [ax S(x— х) (5) Я 
L4 2 / , de E 2U y’ 
and y =V H +30, 4X? — | dx’ 5(х-х) ГЕТ АХ”). (16) 


Thus we get а nonlinear non-local Schróinger equation 


Wee y [ar еа) (a^) = 0. (17) 


In the following we shall first evaluate the correction to the dispersion relation 
(5) due to the nonlinear effect of the amplitude. Then the analysis described 
above will immediately lead to our desired result. 


4. Correction to the linear dispersion relation 


In the present case of a relativistic plasma, we go back to (2). To obtain the 
desired nonlinear correction to (5), we take spatial Fourier transforms of f,, f, 
and ¢ and Laplace transforms in time (Roy Chowdhury & Roy Chowdhury 
1991; Chakravorty, Pakira & Roy Chowdhury 1993): 


f.(K,t,v,) = [е t,v,) e !F* dx 


f(K,S,»,) = | Каз)“ Res >S, (18) 
0 
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the Laplace inversion being 
1 Sotia 
cA. st 
fe (h.t, Ve) Omi Lagar tee 48, (19) 
where it is understood that an appropriate semicircular contour has been 
chosen. Similar formulae also hold for f; and ¢. 


We can now rewrite (19) in а more compact and useful way. We set S = — iw; 
then 
1 tt ia. 
f(x, v, t) =—— do | dK (К, о, о) 82-90, (20) 
(2)? J а 


so, instead of (3), we now substitute 


ee Fe) [fp vg eee do (21) 


where р stands for (К, о) and dp = —dK do/(2n)*. Similarly, we set 
fix, va t = Fav) + [део e Emod dp, 
P(x, t) = fæ et (E ot) o(p). (22) 
Equation (??) yields, 


JAP, Ve) = (w— ko,)! fie K,v,,t = 0) Ky v.) plp) Fy (о 
г, —(K—K')v n v,,t 0) 
+ xt) (К-К) $(p — p’) Fi(v,) 


+= fison K'D,f(p-p' p", 2! (23) 


In the Laplace convolution integral we choose the w contours to lie above all 
singularities of ¢(K,w), and the w—w contour to lie above all singularities of 
f(K,v,w — w'). Iteration of (23) leads to 


fip.) = (p) + o- Kop? [Eyed dtr) 
jore dco A 


eJ [apap KoD tte-en- ica ne 
xD (x) ККК") $(p—p' -p ша в”) 


or (K—K' =K") PET (24) 
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where I (р, v,) consists of all the inhomogeneous terms, namely those linear in 
the initial perturbations, and each operator D, acts at all functions of v, to its 
right. A similar analysis of the ion and Laplace equations yields 


Л, vi) = Hp, v) =w- Kw) ков) 
- арк фр ow = (KK uP K-K) др?) 


2 
(E) [ae ai к агро КК), 
4 


x (К-К – К") (р-р —»')[(»—o—o)—(K —K' —K") DELI 


(25) 
—К*$(р) = —4mem, [ dv, + 4тет, Trento, va dv,. (26) 
With the help of (24)-(26), we get 
—K'e(p) ф(р) = bop) + | dp’ [m М, (р, р’) +m, My, (p.p) Op’) (р-р) 
+ је dp" [m М; (p, р’, p^) 
+m Map P pP NAPA) p-p р"), (27) 
where e(p) = a i-em [im = 
_ Х(о,) (Ve) Fy (Ue) 
47e? 0:5 (бшу) a», (28) 
Фо(р) = т, neat. Ve) do, —m, [м v) до. (29) 


In (27) and (29) the symbols M,,, Mio М», and M,, stand for the following 
complicated expressions involving the initial distribution functions f and x, 7: 


My(p,p) = (=) ккк-к) IMEEM (30a) 
mop (5) кокко [2627 Теа уюду, Bo 
М (р, p.p) == (=) кккчк-к-к” 
le dm Tomy 009 
оо (КК) о] (оо 7 (КК Куо] 
Map. p.p) -(£) eae m dv, ens e| 
D, Gro) Fi (n) "T 


* lo—o)—(K—K)v][(o—o —')-(K -K K)o] 
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Equation (27) is now used to find the nonlinear frequency shift. Since this 
equation cannot be solved exactly, we again make use of perturbation theory. 
Let us consider the iterative solution of (27) in the following form, up to terms 
of third order: 


Wig’) AD (a — 9’ 
$t) =- Jar fm My or т.м р РР, (з) 
АЗ ét) = [ap mM, i) +m М рор SM) doo ph 


+ 17 p)tm,Ma(p.p)9? (p) $9? (р-р) 


+ far dp” [mM (p, p', p^) 


+т„М (р, p,p NAP) 99 (p") 9 (p—p'—p"). (32) 


We are interested in the case where Q(K, w) has a peak value. Let ф(К, w) be 
peaked around К = t К; then we can set 


__i6({K) 
$(K, w) zd w—Q(K) > (33) 


where Im (v — O(K)] > 0. For the above iteration process, we take 


:$(К) 


g® = meray (34) 


It is then not difficult to see that $(2) is zero: if we demand that ¢" and its 
iterated values should all be peaked near К = + К, then, from (31), we find 
that all three modulii |p|, |p’| and |p— p'| should be equal to К, ог —K, 
simultaneously, which is impossible. The next finite contribution comes from 
the third-order iteration ¢@. 

We now choose the w’ and w” contours in (32) to be infinite semicircle in the 
lower half-plane; we then enclose the poles at о =Q(K’)=Q’ and 
ш” = Q(K") = О” coming from the factors of ¢6%(p), and those at w = K'v, of 
Mj, (p, р-р) and о = K'v, of M, (p, р-р). From (32), we immediately get 


— K'e(p) 9p) = Emi акаа", o") [m МККК", 


m, MG, Q0; K’, ; K”, Q”) 

—([m;M,,(K, 0; K —K', Q—Q") -m,M,(K, 0; K — K', Q—Q')] 
—[m,M,,UGQ ; K', ^) -m, MG О; К’, О} 

x [mM (K-K Q—Q' ; К", Q") -m,M,(K —K', Q—Q'; К", ony} 


x (K K')* «(p —p')? $9 (K', Q^) ф®(К”, 0") 49 (K — K' -K",Q—Q' —Q"). 
(35) 
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We now define e(K) = e(K, Q(K)), 
R (K, K, K^") = M,,(K,Q; K,O'; K",0^), 
R,.(K, K, K^) = M,,(K, Q;.K’, О’; K^, О”), (36) 


R,(K, BK’) = My (K, Q; K’, Q’), 
Р, (K, K’) = MAGO, Q; K', О’). 
Using these, (35) can be rewritten more compactly as 


— K'e(p) 99? (p) 
= (2ny?i fax dK” R(K,K', K^) ф(К') $k") (K -K' К”), (37) 


where the kernel R(K, K', K^) has the form 
R(K,K', K”) = m,R,,(K, K', K^) -m, R,,(K, K', К”) 
— [m, R,,CK, K — K^) - m, R,,(K, K — К’) 
xm, RGGKD) +m Ro (K, K)] [m R,,(6 — К’, К”) 
+m, Ra (K-K; K)) (K К) ?[e(K -K)] ! [o Q' - O0" —0(K - K' - K^)! 
tT(K, K', K^). (38) 
The term 7(K, К”, К”) is 
T(K,K', К") = - ККК — K') [m M (K-K, o —K'v,; K",Q") 
+m,M,,(K—K’,w—K’v,; K",Q")] 


E d. dv, Р, 
т, J (o — Ко) (o К'о, —Q" - Q(K — K' - K")] 
x (K'v, —Q' — iy) [e(K -K',o — К)" 
EE 3 , 7(.) К 7 (ve) Ea -O”— 4 К" 
n [ехе | ee 97 (o — K'v, —O" —Q(K — K' — К”)] 
х (K'v,—Q i) ! [(K -K' ;o -K'v)] (К-К), 0,9, > 0 (39) 
This term actually comes from the poles at w = K'v; and о’ = K'v, in M,, and 
M, respectively, but makes no contribution to the final result since it has no 
factor (w—Q)7}. 
We now express К, К”, and К” in terms of small deviations 4, 4’ and q” from 
K,. We write 
Ry (K,K’) = Ra (K, Ko +7) 
= MQG,QO; Ky 9, Q(K, 9), ete. (40) 
Since ¢(K,w) has a peak near К < +К,, we can assume that the initial 
perturbations f,(K, ve = 0) and (К, v, t = 0) are non-zero only for К z t K,. 
We then assume that the wavepacket can be written as 


dz, t) = p(x, t) согоо + С.С, 


with ylz, t) = (27)! [м у (д) e», (41a) 
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with v = Q(K,-- q) — Ә(К,), whence, from (35), we get 


— Кф (p) e(p) = (21) * i IE dq" 


XS q.d) WO) V*(-a") V(q—a —)) (e— Ә(К)Г", 

S(q, q^, d")  R(K, +g, Къ+9, — Ko") E R(Ko rq, Ko tg, Ko +9”) 
+R(K,+9, —Ka +g, Ky q^). (416) 
The explicit form of R can easily be ascertained from (37). The expression is 


very long, and we shall only give here a term representative of its structure. For 
example, 


€ 3 , "n ГА ГА 
В, (К +9,Ко+9;К) --(£) (К, +9) (К,+9) K"(q —q' —q^) 


x одак, догоо) ЧФ", (42) 


with similar expressions for E,,(K,-- q, K,-- q^, K”), R,í(Ko4- q, Къ+9, K^) ete. 
Setting q = q' = qg” = 0 in all the terms without poles, we can at once obtain 
simplified expressions for them. 

There are many terms with w— kv; and w— Kv, in the denominator, which are 
handled by the usual Landau damping procedure. For this, let us consider the 
term R}: 


В (К, K’) = (=) кк(к-к) fanr; [0—0 — (K — K') v] ! (o — Kv) ?, (43) 


lim к-к)|& Мо) ЦК К) (vo — Q')] 


к-к 


> P |а MY) арз (K—K’)ih(v,) (44) 


0 — 00 


with w = О +1, > 0,, and sgn (x) = x/|x|. So we can let К, K' ^ К, in R,,, but 
the result is still a function of K' — K, with 


K=K,+q, К’=К +9’, (9,9,9—9,9') 0. 


Finally, we obtain 
В (К,К)> oy (K,,q—94), 


"NS -(£) dv, F; 
TT Th (v; —9, — 0/9) (v, — u)? 


with 


= (К) +в (9) 0K), (45) 
e F, 
“= E А оао), 
Fi.) 
w = п 
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Similar contributions from other terms lead to 
В. (К, K’) > @(Ko,9-7) 
= (К) +i sgn (д) 04(K), 


os) = (LY Pl (an, F9 x9 2 
т, 


e. d (v, — Vg) (v, — uy? 


n | an Feo Qaa] 


(v, — Vg) (v, — и) 
ev 
w(K) = (2) п [EF (v5) (v5) (о) (v, — u)? 
FQ) 9 (3) Xg) (v, — w) 1), 
Ry(K, Kyo tg, Ky q^) - E (К, Ky td, Kg 07) > Cu (K9,9— 4), 
where С. (K,q) = С, +С, sgn (q), 
е 


С. = -(£ P fan (v, u)? (v, — v) D Р.о —v,) (v, — м), 


m; 
e 3 
С, = -(GJ nF (v) (v, —u) *, 
Ry (K, Ky 7, Ky 9") + RC Kyo tq, —Ky tq”) > Coo(Ky, 9-7), 
C..(K, q) = С, +С, sgn (9), 


C, = EO, PÍ fnt) (v, — 4) (9,9) " 


€ 


x D, [x(v, — Ug) Fy (ve =0,)] (v, pm u)? 
F fev. 7 (v, xS Ug) (v, Re u) (v, Ж Oye 
Dit) Во, 910—307), 


eY Ж _ а 
C, = -(J п [x(v,) F1(v,) (v, — и) ^ xlv) F,(v,) (v, — и) =] . 
Furthermore, 

Ry (K—K’, Kor q^) +В (К-К’, - Ko 9") > wn Kogi), 
where w,, is given by (45). We also observe that 


R,(K—K’,K,+q')+Ry(K—-K’, —Ky+9")> 0.:(Ky, 9-7) 


(47) 


(48) 


(49) 


(50) 


where w,, is given by (46). A contribution to the singularity in К — K' also comes 


from (K—K’)* є(К — К’). We have lim (K—K’)? e(K—K’) = —^(K,9—4), 
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where A(K,q) = [4 on -Ay + [aries —% 28). 
= A,(K) t isgn (9) A(K), 
A, (K) = P| | а,в нз |а, ey (51) 
АК) = по) - Fit. (52) 


5. Contribution of non-singular terms 
The contributions from terms without any singularities can easily be 


evaluated by setting = = 0 and К = К, in the expression (416) for 
S(q,q', q"). We have 
3 
R,Q(K,K', K^) = -(&) KK'K"(K — K' К^) fos (o — Kv)? 
4. 
х [0—0 —(K —K')u]? 
хр, Е (и) [(» — »' —^) — (&—K' K”) v], (53) 
31 
R,,(K, -K, K) = (=) |е (v, — 4) *D, Fy(v;) (ии), (54) 
i 
=4B_ (say). 
Similarly, from (304), 
вык, К.К) = (Z) 5} [advent бе Ditta Fea и 
x duet (ve) (e, — ро Fi) Weu} = 0 (вау), (55) 


Е (К, -K) = (=) [ario (—w)? = -А (вау), 
Р, (К, —K) x (еледе X(v,), (v, —u)? 


zd Jo. 3j (Ve) F'i(v,) x(v,) (v, Ze 04 


= -Е (say). (56) 
Other relevant contributions are 
RQ(K,2K)-2A, R,,(K,2K) = 
R,,(2K,K)=14, R,,(2K,K) = 15, 
R,,(K,K-K)=0, R,,(K,K—K) =0, 
T(K,K, +K) = 
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(57) 
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Substituting these into the expressions (416) for S(g,q',q"), we get 


S(q,9',9°) = S (K) —isgn (9—9) 5(К), (58) 
2 
S, АВЕ — Wm Bem, D) m, C, tm, C, 


+ (о, о +m, ws)? — (m, v, +m, в. ^, 


t 2(m, ө, t m, Ws) (Mi w + m, w4) Ag} (AF + A3), (59) 
S, = — (m; C, m, C4) — (2(m, ш, +m, ws) (m, v, +m, v4) ^, 
— [(m, оу +m, ws)? — (m, w +m, о,)?]). (60) 


6. Derivation of the non-local nonlinear equation 


We now proceed to derive the nonlinear equation after the corrections to the 
linear dispersion relation have been evaluated as above. Let us now consider the 
basic integral equation (27) and recall our basic assumption that the initial 
perturbation f(K,v,t = 0) is non-zero only for К »+K,, that is, when ¢ is 
peaked around K,. To make this explicit, we set 


К) = (К-К) +*(—K—K°), (61) 
with (9) non-zero only for 4 z 0. Then (27) becomes 


—K'e(p) $(p) = Фр) + i(27) [жа 


x Sig) y(a—a') 009 —4)у*(—4')[в—О(К)|1, (62) 
where we have changed the integration variable from q' to q—q' so that 
S(7) = $(0,4—4,0) = S,(K) — iS, (K) sgn (q^). (63) 
Putting the same form of ф(р) on both sides of (62), that is, 
$(p) = К-К) c y*( -k- K,)] lo- QK)”, 


and considering the region of K near K,, we arrive at 


(p) $(g) + (27К)- fay dq" 


харуу aov) ilo) fh. (64) 


If the second term of this equation is proportional to (4) then the coefficient 
of (д) will give the frequency shift required for (17). This proportionality is 
most easily displayed through the Fourier transform with respect to q by 
treating p — [K, O(K)] as constant. For the left-hand side of (62), we then obtain 


[en [зрение e] 
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which has the same form as (17). Therefore it is now possible to identify the 
kernel Q(x) in the NLS equation: 


S, [1 
дш) = | 0e) +p (2), (65) 
whence the non-local nonlinear Schródinger equation becomes 
d 
iht phata ype ур [MEN _ 6, (66) 
_ 8, (0e? _ 8, (ôe)! 
where Q = az) ‚ Ф = (с) ; (67) 
де (о) dv; | ((v,) n(v,) Е! (в) dv | 
2| "" 21 e e 1V^e [4 
К (=)=- «el [rum UE v,—w/K ` (99) 


7. The case of a Maxwellian plasma 


Although it is known that the Maxwellian distribution does not apply in the 
relativistic situation, we shall still use it for the explicit estimation of the 
parameters in (66), since none of the integrals can be obtained exactly with a 
modified distribution function. Let us assume that the unperturbed velocity 
distributions for the electrons and ions are 


Бы) = тер -5(55) | 
(69) 
Fu) = nmiexp| -5(5] | 


where the electron thermal velocity 
and the ion thermal velocity 


The relation between the frequency w and wavenumber K turns out to be 


3K? xs 
w? = Mm, We - m, wp? + pe (e Оре Vo m Opi Va) 


4 


K 
+ 15-1 (Me Ope Vo +m, wiva) (70) 


4me^n 
where Upe = ——, ар; = 
m, т 
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are respectively the electron and ion plasma frequencies. Now we can evaluate 
the integrals with the help of (69), and we obtain 


2 К? 2 1 4 
w,(K) = -(5 Arm и (71a) 
m, wv? vi v; 
е\ К? о, 
w,(K) = (=) LO (715) 
t 
2 2 2 
w(K) = (5) Е K piae 
т, vw v v 


ЗК [за A ть 
+ Co s(t % 


К 
n) (5+ СИ d (114) 


(71с) 


т,) Vo? 
‚=-(& Ja lat qs oe em), (72a) 
e К, 
С, = -(£) = — Felt) (726) 
eV ГК ЗК 15K*wj 
&--() | art oe ETE 
е, 9 g g 
КЗ (25$ , 5%\ 6K? 
seal vi % v? tons ў 220 
аъ КА, 
С, = -(J uri (72d) 
1 2 1504 1 2 15v; 
^ 2 ек І). (73а) 
Vg Vg 9 о а % 
A, = п([Е,(о,) +Р.(0,)]. (735) 


Finally, v, = ды/0К can easily be computed from (70). With these expressions 
at hand, we now have explicit forms of the coefficients Q,, @„ and р occurring 
in (66). These are rather long, and so are given in the appendix. 


8. Nonlinear Landau damping 

We now proceed to estimate the amount and variation of nonlinear Landau 
damping by an approximate analysis of (66). Let us consider a perturbed 
solution of (66) in the form 


where p and c are real functions of £ and 7. We set 
р = р+р; ei(KE-0n). | (75) 


о = og, o, 606-97), ] 
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FiGURE 1. Wave vector given by the dispersion relation as a function of the frequency 
ratio о/о р. for v, = 0-5, w = 0-5 and different vp. 
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Етасве 2. Wave vector given by the dispersion relation as a function of the frequency 
ratio v, for w,,/w,, = 2; ө = 0-5 and different v,. 


where p, < p, and c, < с,. Substituting into (66), we get 


(Q—o, K) = — K?p (2p, [Q; —i sgn (K) Q, — [р]). (76) 
Setting О = AW cil, Юг К > 0 we get from (76) 
T = + Gp (20, Qi — DK?) + (20, Q} E+ (20, Q, — pK?) К. (77) 


which gives the nonlinear Landau damping as a function of several plasma 
parameters. To ascertain the behaviour of Г, we start by analysing the 
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Ficure 3. Nonlinear Landau-damping rate Г as a function of frequency ratio w,,/w,, for 
v,/¢ = 0.5, w = 0-8 and different v,, exhibiting the tendency of growth for small v, and large 
w. 
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Ficure 4. Nonlinear Landau-damping rate Г as a function of frequency ratio w,,/w,, for 
v, = 0-5, w = 0:5 and different v,/c, exhibiting the tendency of growth for large v,/c and large 
Ф. 


dispersion relation (70) for several values of v,, as a function of w,,/w,, and v,. 
The corresponding solutions of the wave vector K are shown in figures 1 and 2. 
Equation (70) is then used to estimate the behaviour of the nonlinear Landau- 
damping coefficient Г. This is clearly shown in figures 3-6. In figure 3 we observe 
that the Landau damping coefficient Г grows quite steeply with w,,/w,, for fixed 
v/c and о but different v,. In figure 4 we show the variation of Г for fixed v, and 
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Боске 5. Nonlinear Landau-damping rate Г as a function of frequency ratio w,,/w,, for 
v, = 0:5, о = 05 and different v,/c, exhibiting the tendency of growth for large v/c. 
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Ficure 6. Nonlinear Landau Damping rate Г as a function of the ion thermal velocity v, 
for different v,/c, but fixed w,,/w,,. 


w but different v,/c. For low values of v/c the damping coefficient reaches a 
constant value, but for larger v,/c it again increases sharply. On the other hand, 
in figure 5 the same variation is shown for lower frequencies. The qualitative 
behaviour remains the same, but the magnitude of Г is scaled down by a factor 
of 20. In figure 6 the variation of T with respect to v, is shown for fixed w,,/@5- 
and о but different v,/c. It is found that the larger the value of v,/c, the lower 
is the value of Г, but ultimately for large v, they converge to same value, almost 
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FrevnE 7. Nonlinear Landau-damping rate Г as a function of v,/c for w,,/w,, = 2, w = 0:5 
and different v,. 


zero. Finally, in figure 7 we have interchanged the role of v, and v,/c; that is, 
we have plotted Г as a function of v,/c for different v,. The same phenomenas 
are seen to occur. 

Thus, from the non-local nonlinear Schródinger equation, we have obtained 
an estimate for the nonlinear Landau-damping coefficient for several values of 
the plasma parameters v,/c, v,, Wpi/Wpe and w. Whether such an equation can 
sustain any soliton-like solution is yet to be seen. 


9. Discussion 


In the above analysis we have shown how the Vlasov plasma kinetic 
equations can be effectively used in conjunction with the approach of Fried & 
Ichikawa (1973) to derive a generalized non-local nonlinear Schrödinger 
equation. It may be noted that a similar equation has been derived previously 
for a non-relativistic plasma by a completely different method. In this respect, 
it may be mentioned that the actual velocity distribution for the relativistic 
case is not a Maxwellian one but rather is a Juttner function. However, none 
of the above integrals could be performed explicitly for such a distribution, so 
we have used the Maxwellian function to describe the plasma approximately. 
Another approach to the derivation of (66) utilizes a Fourier-mode expansion, 
which is more complicated from a computational point of view. The last part 
of our analysis has involved explicit computation of the nonlinear Landau- 
damping coefficient for various plasma parameters, all of which have been 
depicted graphically. 


One of the authors (J.M.) is grateful to the University Grants Commission 
Junior Research Fellowship in a research project. 
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Appendix 
The explicit forms for Q,, Q, and p in (66) are as follows: 


S, (0e! 
а-а) 
1 
= mm Cy +m, C3 + {[( v +m, w3)? — (m, v; + m, 4)"] A, 
t 2(m, оу +m, Ws) (m, ® + m, w4) Ag} (А? + х] 


2 K2 4774 
s EZ (+ ee 


pe 3 "ma "P 


1 602 К 45K" 
+m (+ 1+ ) ; 


wW w 
| S (0e)? 
Q - e) 
1 
= zl- (m, C, +т, C4) — (2(m, 0, +m, оз) (m, Wg + m, W4) ^, 


— [(m, w; +m, w3)? (то, +m, «,)*]) 


2 K2 4 Kga 
x [2m vt РУ я 


e PEN Gy 5 ө 


1 62K? 450% KAT! 
+ 2m, о (5 + wo + ИШ ) , 
is 1д%ф 
Р = 22° 
= B 
~ 9[605 — 4(т, о, +m, о?) 03) 
where 
B = [60° — 4(m, v2, +m, v) о] 
x (ЖКХ РУ | 6x? + 12Ко+ (12kw+ 6K? ac 
+ 180K?(m, 05, 00 +m; о v; 
— [(6kw? + 6K?w) (m, о, ve +m, v2, v?) + 60K? (m, о, vo +m, Wh; V9) | 
x [3004 — 12% (т, о, +m, v?,.)] ш 
ере £^ pi OK" 
with 
Oo _ 6Кш(ш+ К) (m, o5, vg + m, v5, Vz) + 60K? (m, Ope v + m, 95, v;) 
aK бо — 4(m, о, +m, о) e? | 
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